Fragmentation of electric dipole strength in = 82 isotones 
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Fragmentation of the dipole strength in the A'' = 82 isotones ^^"Ce, ^*^Nd and ^*^Sm is calcu- 
lated using the second random-phase approximation (SRPA). In comparison with the result of the 
random- phase approximation (RPA), the SRPA provides the additional damping of the giant dipole 
resonance and the redistribution of the low-energy dipole strength. Properties of the low-energy 
dipole states are significantly changed by the coupling to two-particle-two-hole (2p2ft) states, which 
are also sensitive to the correlation among the 2p2h states. Comparison with available experimental 
data shows a reasonable agreement for the low-energy El strength distribution. 
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The low-energy dipole states, often referred to as the 
pygmy dipole resonance (PDR) , have attracted recent ex- 
perimental [H-Q s-nd theoretical interests (see also 
the recent review [l^ and references therein). It is also 
of significant astrophysical interest, since the low-energy 
dipole strengths close to the neutron threshold strongly 
affect the astrophysical r-process nucleosynthesis [ij] . 

The quasiparticle random-phase approximation 
based on the Hartree-Fock-Bogoliubov ground state 
(HFB-fQRPA) has been extensively used to study the 
PDR as well as the giant dipole resonances (GDR). 
Recent systematic calculations Tsj for the Nd and Sm 
isotopes show that although the HFB-I-QRPA nicely 
reproduces characteristic features of the shape phase 
transition in the GDR, it fails to produce the low-energy 
dipole strengths at iJj; = 5.5 ^ 8 MeV, observed in 
the N = 82 isotones, ^^^^d and ^^''Sm [U, [|. The 
disagreement suggests that the coupling to complex 
configurations, such as multi-particle-multi-hole states, 
are required to study the PDR in these nuclei. In fact, 
the quasiparticle-phonon model (QPM), which takes into 
account coupling to multi-phonon states, successfully re- 
produces the low-energy dipole strengths in the TV = 82 
nuclei 0, 13] ■ A similar approach based on the relativistic 
mean-field model has also been used to study the PDR 
in the tin and nickel isotopes These models assume 
the multi-phonon characters of the complex states and 
violate the Pauli principle. Thus, it is desirable to study 
properties of the PDR with a method complementary 
to these phonon-coupling approaches. In this work, we 
present studies for the dipole excitations in the TV = 82 
isotones, with the second random-phase approximation 
(SRPA) (Ref. [13 and references therein). The SRPA 
explicitly incorporates the two-particle-two-hole {2p2h) 
states instead of "two-phonon" states, and respects the 
Pauli principle in the 2p2h configurations. Recently, the 
low-energy dipole states in ^°''**Ca have been studied 
with the SRPA [l^, which suggests that the couphng 
between one-particle-one-hole {Iplh) and 2p2h config- 
urations enhances the electric dipole {El) strength in 
the energy range from 5 to 10 MeV. We investigate 
whether a similar effect can be observed in the isotones 



of = 82. Since there are many dipole states with 
small El strengths in the energy region below 8 MeV, 
it is difficult to compare property of each state with 
the experiment. Thus, we perform the comparison of 
integrated properties at low energies. 

The SRPA equation is written in the matrix form[l3] 



a c 
b d 



(1) 



where and X, 



pp'hh' 



{p O h) are the Iplh and 2p2h 
transition amplitudes for an excited state with an excita- 
tion energy w^. The explicit expression for the matrices 
a, 6, c, and d are given in Ref. [19]. 

The Skyrme interaction of the SIII parameter set is 
used to calculate the Hartree-Fock single-particle states. 
The continuum states are discretized by confining the 
single-particle wave functions in a sphere of radius of 20 
fm. Single-particle states with the angular momenta ja < 
15/2 up to 30 MeV in energy (ea < 30 MeV) are adopted 
for the Iplh space (a;^^ and x'^^), both for protons and 
neutrons. This roughly amounts to one hundred single- 
particle states. For the 2p2h ampfitudcs (Ar^^,^^, and 
■^hh'pp')^ we truncate the space into the one made of 
the single-particle states near the Fermi level, the 2p2/2, 
2pi/2, 159/2, 157/2, 2(^5/2, 2^3/2, 3si/2, l/iii/2, and 1/19/2 
orbits for protons and the 2d^/2, 2c?3/2, I/111/2, 1^9/2, 
2/7/2, and 1^13/2 orbits for neutrons. The proton orbits 
up to the I57/2 orbit are assumed to be fully occupied in 
the ground state of -^^"Ce, while the proton 2^5/2 orbit 
is to be partially occupied in the ground states of ^^^Nd 
and ^^"^Sm. The numbers of Iplh and 2p2h amplitudes 
in the SRPA are about 800 and 9000, respectively. 

For calculation of the SRPA matrix elements, we em- 
ploy a residual interaction of the and terms of the 
SIII interaction. Since the residual interaction is not 
fully consistent with the one used in the calculation of 
the single-particle states, it is necessary to adjust the 
strength of the residual interaction so that the spurious 
mode corresponding to the center-of-mass (COM) motion 
comes at zero excitation energy in the RPA. This condi- 
tion determines the renormalization factor / for the resid- 
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FIG. 1. Strength functions calculated in the SRPA (solid line) 
and RPA (dotted hue) for ""Ce. An artificial width T = 0.5 
MeV is used for smoothing. See text for details. 



ual interaction {to ^ f x <o and —i' f x t^) . We obtain 
/ = 0.73 for ^"^^Nd, and similar values for other nuclei 
as well. Since the coupling between the spurious COM 
motion and 2p2h configurations is weak, these renormal- 
ization factors may approximately produce zero energy 
in the SRPA as well. Thus, we use this interaction for 
the calculation of the matrices a, b, and c in Eq. ([IJ. 
For the residual interaction for the matrix d, following 
a prescription in Ref. [l9j . we introduce a zero-range in- 
teraction vod^{r — r') in addition to the original and 
is terms, then, fix the parameter vq by approximately 
reproducing the excitation energy of the lowest 1^ state 
in i^'^p^d ^ „57o MeV fm^). With these residual in- 
teractions in the given model space, the spurious mode 
appears at a small imaginary energy (w^ w — 1 MeV^) in 
the SRPA. 

We first show the results for the GDR. The El 
strength functions, S{E) = J2„ \ {n\rYi^\0)\^6{E - E,,) ^ 
dB{El;l- -> 0+)/dE, calculated in the SRPA (solid 
line) and RPA (dotted line) for ^^°Ce, "^Nd, and ^-^-^Sm 
are shown in Figs. [Tl[5]and|31 respectively. We use the El 
operator with the recoil charges, Ne/A for protons and 
—Ze/A for neutrons, for the calculation of S{E). The 
obtained discrete strength functions are smoothed with 
a small width (F = 0.5 MeV) of Lorentzian. The energy- 
weighted strength summed up to 50 MeV exhausts 87% 
of the energy-weighted sum-rule value including the en- 
hancement term arising from the momentum-dependent 
parts of the Skyrme interaction. The strength distribu- 
tions of the GDR in the SRPA are broadened, compared 
to the RPA, due to the coupling to the 2p2h states. In 
the inset of Fig. [21 the total photoabsorption cross sec- 
tion (solid line) calculated in the SRPA is compared with 
the experimental data [20|. The shape of the GDR de- 
pends on the parameter /, whereas it is little affected 
by the parameter vq- The GDR peak position and the 
profile are better described by a slightly larger value of 
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FIG. 2. Same as Fig. [T]but for ^"^Nd. In the inset the to- 
tal photoabsorption cross section calculated from the strength 
function in the SRPA (solid line) is compared with the exper- 
imental data (dots) [20| . The dotted line in the inset denotes 
the result calculated with / = 0.9. 
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FIG. 3. Same as Fig. [T]but for "''Sm. 

/ (See the dotted line in the inset of Fig. [2]). Our calcu- 
lation indicates that the coupling to the 2p2h induces an 
additional broadening due to the spreading width, how- 
ever, the peak position is close to that obtained in the 
RPA calculation. This is very different from the recent 
SRPA calculation for ^^O in Ref. [23|, which indicates a 
large shift of the GDR peak energy (more than 5 MeV) 
but almost no broadening. At present, we do not fully 
understand the origin of this discrepancy. More quan- 
titative analysis of the GDR require an improvement of 
the present calculation, especially, a self-consistent treat- 
ment of the residual interaction and the enlargement of 
the 2p2h space. 

Next, let us discuss the low-energy El strengths. In 
contrast to the GDR at high energy, the truncation of 
the 2p2h configurations is supposed to be less serious. 
The El strengths, B{E1) i, below 10 MeV in "OCe, 



3 




5 6 7 

E [MeV] 



FIG. 4. Low-energy El strength distributions, B{E1;1 — > 
0+ ) calculated in the SRPA (solid hne) and RPA (dotted line) 
for "»Ce. 
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FIG. 5. Same as Fig. Hbut for "^Nd. 
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TABLE 1. Mean energies E and summed B{E1) t values for 
the low-energy dipole states. The experimental values are 
taken from Ref. [3||. See text for details. 





E^ [MeV] 


EB(£l)t [e^fm^] 


Nucleus 


RPA SRPA Exp 


RPA SRPA Exp 


^^"Ce 
Nd 
"^Sm 


7.53 6.47 6.28 
6.31 6.07 
6.04 5.69 


0.021 0.219 0.308 
0.0 0.224 0.184 
0.0 0.233 0.208 



FIG. 6. Same as Fig. Hbut for ^''^'Sm. 



"'^''^Nd, and ^^'*Sm are shown in Figs. |31[5l andlHl respec- 
tively. In the RPA calculation, there is very little El 
strength in the energy region below 8 MeV, which agrees 
with the result of the QRPA calculation ^5|]. However, 
this is different from the experimental findings [H-ll]- 
the SRPA calculation, the coupling to the 2p2h config- 
urations leads to a considerable El strength in this en- 
ergy region. To make a quantitative comparison with ex- 
periment, the mean excitation energies and the summed 
B{E1) t values for low-energy dipole states are calcu- 
lated in the same way as the experiment The mean 
energy is defined a,s E = EB{E1)/Y^ B{E1), in which 
the summation is performed for the dipole states below 
7.7 MeV for ^^°Ce, those below 7.1 MeV for i^snj, and 
below 7.0 MeV for ^^-^Sm. The lowest Ij" states are ex- 
cluded in the summation. The result is tabulated in Ta- 
blelD For comparison, the RPA values, which include the 
lowest 1~ state, are listed in the table, but no 1^ state is 
predicted below 7.1 MeV for "z^d and "^Sm. Although 
the calculated mean energies are slightly larger than the 
observed values, their isotone dependence is consistent 
with the experiment and the summed transition proba- 
bilities are comparable to the experimental values (Sj] . 

In the RPA calculation, the neutron excitations are 
dominant in the low-lying states 13|. The present RPA 
calculation also indicates, for instance in -"^^^Nd, that 
the largest components of the low-lying dipole states 
located at E^ = 7.36, 8.64, 9.15, and 9.55 MeV are 
(2pi/2 2d3/2)7r, (3si/2 3^3/2)1^, (3si/2 3^3/2)1/, 
and (3si/2 3pi/2)j^, respectively. In the SRPA, we see 
a significant fragmentation of the dipole strength into 
the energy range of 5 < i? < 8 MeV, in addition to the 
emergence of the lowest lj~ state at E « 3.5 MeV. Many 
of these low-lying dipole states consist of proton 2p2h 
characters, such as ([137/22^5/2]^ — > [1^11/22^3/2]'^ )7r 

and ([1.97/22c?5/2]''^ ^ [l/iii/23si/2]''' )'^- These proton 
2p2h configurations come down to the lower energy be- 
cause of the coupling to the 2p2h configurations consist- 
ing of the neutron Iplh transition from the I/111/2 or- 
bit to the lft.9/2 orbit and the proton Iplh transitions 
from the l(?7/2 orbit (or 2^5/2 orbit) to the lft.11/2 orbit; 
■nlg'j/2vlhiii2 -> 7rl/iii/2J/l/i9/2 and. 7r2d5/2J^l/iii/2 ~^ 
"'1^11/2^^1^9/2 and. We have confirmed the importance 
of these proton-neutron 2p2h configurations by perform- 
ing the SRPA calculation in a smaller 2p2h space. The 
SRPA calculation with the neutron Ih orbits qualita- 
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TABLE II. Excitation energies E and B(£1;0+ 1") of 
the lowest 1~ states. Note that the energy of the 1^ state in 
'^^^Nd was approximately fitted by adjusting the parameter 
vo- The experimental values are taken from Ref. 





E [MeV] 


B{E1) t [e^ 


fm"] 


Nucleus 


RPA 


SRPA 


Exp 


RPA SRPA 


Exp 


^*"Ce 


7.53 


3.54 


3.644 


0.021 0.076 


0.0217 


"2Nd 


7.36 


3.45 


3.424 


0.010 0.074 


0.0211 


"*Sm 


7.18 


3.35 


3.226 


0.004 0.068 


0.0248 



lively produces the same result. 

Finally, let us discuss the property of the lowest 1~ 
state. The excitation energies and the reduced transi- 
tion probabilities B{E1) '[ of the Ij" states in ^^°Ce, 
"'^^^Nd and ^^*Sm are compared with the experimental 
values [I in Table HH The calculated excitation energies 
decrease with increasing proton number, which is con- 
sistent with the experiment. However, the SRPA calcu- 
lations overestimate the B{E1) t values by a factor of 
2.7 — 3.5. The structure of the lj~ states in these nu- 
clei is supposed to be predominantly of the two-phonon 
quadrupole-octupole character 2+ (g) 3~ 0, Q- However, 
in the present SRPA calculation, the 2p2h configuration 

{[TTlgr /2vlhii^2V ~^ ['^^hii/2vlhg/2]'^ ) is dominant 
in these lj~ states, which differs from the two-phonon 
2"*" (g) 3~ character. The pairing correlation, which is not 



taken into account in the present calculation, may play an 
important role for a better description of the two-phonon 
character of the l|j~ states, because they are essential in 
the description of the lowest quadrupole and octupole 
states. Furthermore, it has been known that the SRPA 
fails to describe the collectivity of the two-phonon states 
[2^ 1. This is due to the fact that the next-leading terms 
in the two-phonon state are missing in the SRPA. These 
missing terms beyond the SRPA can be taken into ac- 
count by introducing Xphpi^i amplitudes in Eq. ([T]). A 
general equation for the extended RPA formalism with 
the ground-state correlation is given in Ref. [IJ. An- 
other possible method to improve the description of the 
two-phonon states is the dressed-four-quasi-particle ap- 
proach proposed in Ref. [2^1 . These are beyond the scope 
of the present work, but of significant interest in future. 

In summary, the fragmentation of the dipole strength 
in the iV = 82 isotones, "^Nd, "z^d and "4Sm, was 
studied using the second random-phase approximation 
(SRPA). The SRPA successfully produces the spreading 
of the giant dipole resonance and the concentration of 
the dipole strength in the low-energy region, simultane- 
ously. However, the transition strength of the first dipole 
state was overestimated in the SRPA, indicating the ne- 
cessity of a more elaborate treatment for the states with 
the two-phonon character. The calculation based on the 
extended RPA with the ground-state correlations is of 
great interest and currently under progress. 
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